Logarithmically Slow Expansion of Hot Bubbles in Gases 
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We report logarithmically slow expansion of hot bubbles in gases in the process of cooling. A 
model problem first solved, when the temperature has compact support. Then temperature profile 
decaying exponentially at large distances is considered. The periphery of the bubble is shown to 
remain essentially static ("glassy") in the process of cooling until it is taken over by a logarithmically 
slowly expanding "core" . An analytical solution to the problem is obtained by matched asymptotic 
expansion. This problem gives an example of how logarithmic corrections enter dynamic scaling. 
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Dynamic scaling behavior of extended nonlinear sys- 
tems out of equilibrium have attracted much attention 
in different areas of physics II]. In continuum mod- 
els dynamic scaling is intimately related to self-similar 
asymptotics of nonlinear partial differential equations [^) . 
Sometimes logarithmic corrections enter dynamic scal- 
ing laws [jlj. No general scenario for their appearance is 
known. One can say from experience that they appear 
in marginal cases, dividing (in an appropriate functional 
space) regimes with qualitatively different behavior. The 
aim of this work is to investigate one particular setting, 
of general interest, where logarithmic corrections to scal- 
ing appear somewhat unexpectedly: cooling dynamics of 
hot bubbles in gases. 

Heat transfer in gases, strongly heated locally, looks 
quite differently from the simple picture provided by the 
linear heat diffusion equation. The difference is mainly 
due to the small-Mach-number gas flow that develops 
(even at zero gravity) owing to intrinsic pressure gradi- 
ents. This conductive cooling flow (CCF) brings in cold 
gas from the periphery and strongly modifies the cooling 
dynamics. Some aspects of CCFs have been studied ex- 
perimentally and theoretically, mainly in the context of 
the late stage of strong explosions H-p) . In this Letter we 
report a new, striking feature of a CCF. We find that, if 
the initial temperature profile decays rapidly enough at 
large distances [like exp(— fc|a;|) , k > 0], the hot bubble, 
while cooling down significantly, expands logarithmically 
slowly. 

Starting from the continuity, momentum and energy 
equations for an inviscous ideal gas at zero gravity, and 
employing the small-Mach-number expansion, one ar- 
rives at the following nonlinear equation for the scaled 
gas temperature || : 
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where indices t and x stand for partial derivatives (a 
slab geometry is assumed), and v is the exponent in the 
power-law temperature dependence of the heat conduc- 
tivity of the gas M . 



The scaled gas pressure stays constant (and equal to 
unity) in this approximation, so the scaled gas density 
is simply p(x,t) — T^ 1 {x,t), while the gas velocity is 
v(x,t) = T u d x T [|. Therefore, once solving Eq. (0) for 
the temperature, one can easily find all other variables. 

Eq. (hi) has a multitude of similarity solutions: 



(x,t) = t^6(x/t^), 



(2) 



where (3 is an arbitrary parameter. Therefore, an inter- 
esting selection problem appears, like in many other sit- 
uations in nonlinear dynamics of extended systems . 

Eq. (|]) has appeared in the context of cooling of the 
"fireball" produced by a strong local explosion in a gas 
An explosion involves energy release on a time 
scale short compared to the characteristic acoustic time. 
In this case the preceding rapid stage of the dynamics 
produces an inverse power-law dependence of the gas 
temperature on the distance from the explosion site [||. 
It has been shown that the exponent of this power 
law uniquely selects the scaling exponent [3. As the re- 
sult, the fireball expansion exhibits a power law in time. 

A different type of local gas heating occurs when the 
time scale of the energy release is long compared to the 
acoustic time, but still short compared to the cooling 
time. In this case the initial temperature profile is more 
localized, as it reflects the spatial structure of the heating 
agent (for example, the radial intensity of laser beam). 
This regime will be in the focus of this Letter. We will see 
that there is no self-similar asymptotics to this problem. 
Instead, the solution approaches, at long times, a "quasi- 
similarity" asymptotics with logarithmic corrections to 
scaling. 

Consider first a model problem when the initial tem- 
perature profile has compact support: T(x, 0) = Tq(x) > 
at x E [-L, L], and zero elsewhere. We will limit our- 
selves to a temperature-independent heat conductivity, 
v = 0. Despite this choice, the nonlinearity of Eq. (|f|) 
persists. Assume symmetry with respect to x — and 
impose the Neumann boundary conditions: <9 x T(0,i) = 



f 



d x T(L, t) = 0. A local analysis of Eq. ([!]) near the edge of 
support of its solution T(x, t) shows that the support re- 
mains compact and unchanged for t > 0. Therefore, this 
model problem exhibits a complete localization. What 
is the late-time behavior of the temperature? The con- 
stancy of support immediately selects j3 — 0, so the sim- 
ilarity ansatz becomes To(x,t) = t~ 1 6(x). Then Eq. (|l|) 
yields 6{x) = (a 2 / 2) cos 2 {x/ a) for x G [-L,L],6(x) = 
elsewhere, and a = ttL/2. This simple similarity solution 
describes cooling of the hot bubble (and filling it with 
the dense gas) without any change in the bubble size. 

Remarkably, To(x,t) represents a long-time asymp- 
totics for any initial condition that has compact support 
[—L,L] and obeys the Neumann boundary conditions. 
We will show here only that this solution is linearly stable 
with respect to small perturbations, and find the spec- 
trum of the linearized problem. Introduce new variables: 
u = t,T(x, t) and r = Int. Eq. (|l|) assumes the form 

d T u = u — (d x u) 2 + ud xx u , (3) 

while the similarity solution To(x,t) becomes steady- 
state solution 9{x). Introducing a small correction v(x, r) 
to this solution and linearizing Eq. (||), we obtain 
d T v = where 

L(f) = (1/2) cos 2 £ % + sin 2£<% + 2 sin 2 f (4) 

and £ = x/a. We look for the eigenfunctions in the form 
of v(£, t) — e 7T, 07(C)- The general solution of the result- 
ing ordinary differential equation is 

t/> 7 (£) = Ci cos 2 C2i ? i(a-,a+,l/2,-tan 2 
+ C 2 cos £ sin £ a-Fi (6- , b+, 3/2, - tan 2 , (5) 

where 2 -Pi is the hypergeometric function, C\ and C2 
are arbitrary constants, a± = [l ± (87 + 9) 1 / 2 ] /4 and 
b± = [3 ± (87 + 9) 1/2 ] /4. Requiring that the perturba- 
tion remains small compared to the unperturbed solution 
[and hence vanishes like (ir/2 — £) 2 or faster at £ — + tt/2], 
we find the (continuous) spectrum of the linearized prob- 
lem: — 00 < 7 < — 1. This result proves linear stability 
of the similarity solution T (x, t). Notice the presence of 
gap between the upper edge of the spectrum 7 = — 1 and 
stability border 7 = 0. Going back to physical variables, 
we find that small temperature perturbations around the 
similarity solution exhibit a power law decay i 7_1 . 

At this stage we notice that (3 — is a marginal case 
dividing two qualitatively different types of dynamics as 
described by the family of similarity solutions (||). In- 
deed, solutions with (3 > describe power-law expansions 
, while solutions with p < correspond to power-law 
shrinkings [Qj . One can expect logarithmic corrections to 
appear "on the background" of the special case /? = 0, 
when the initial condition does not have compact sup- 
port, but decays rapidly at large distances. Therefore, we 
assume that the initial temperature profile of the bubble 



is symmetric with respect to x = and decays expo- 
nentially at large distances (To). We will continue using 
the new variables and require u(x,0) — > cexp(— k\x\) at 
|z| — > 00, where k and c are positive constants. One 
can always put k = 1 JO]. We will be interested in a 
long-time asymptotics of the solution: t> 1. Our first 
important observation is that u(x,t) — cexp(r — x) is an 
exact solution of Eq. (^|). This traveling wave solution 
with a unit speed corresponds to a steady-state solution 
T(x) = cexp(— x) in physical variables, and it represents 
the correct asymptotics of the solution to our problem at 
x — > +00. What about the bubble "core"? We will show 
that it can be described, atr> 1, by a "quasi-similarity" 
solution plus small corrections: 

u(x, t) = u (x, t) + Ui(x, t) + . . . , (6) 

where 

uo(x,t) = cos 2 (7) 

2 a(T) 

and ... -C Mi <C u Q . One of our goals is to find an 
asymptotic expansion for a(r). 

The leading term of a(r) can be guessed immediately. 
Indeed, expansion of uq(x,t) in powers of x — 7ra(r)/2 
near the point x — 7ra(r)/2 begins with the term 
(l/2)(x-ira(T)/2). This is a wave traveling with speed 
7rd/2 along the x-axis! Therefore, it is natural to look for 
a general traveling wave solution v(x, r) = V{x — t) of 
Eq. d|) with a unit speed and require that it behaves 
like (z + const) 2 /2 at z — > —00 and like c exp(— z) at 
z — > +00, where z — x — t. If such a solution exists, we 
can match it with the leading term of the quasi-similarity 
solution (Q) in the region 1 -C — [x — na/2) «a,l« — z, 
once 

a(r)=2T/7r. (8) 

Eqs. (0) and (||) have important implications. First, 
the temperature scaling with physical time t at the bub- 
ble center acquires a logarithmic correction. Second, the 
bubble core expands logarithmically slowly. We will show 
in the following that these are indeed correct results, cal- 
culate the subleading and sub-subleading terms for a(r), 
and find other attributes of asymptotic solution. 

The general traveling wave solution of Eq. (||), V(z), 
obeys the second-order equation 

-V Z = V-V? + VV ZZ (9) 

that is soluble analytically. One integration yields 

V- l {dV/dz) = -1 - W (-exp(-l - V- 1 )) , (10) 

where W(rf) is the product log function defined as the 
solution of equation We w = r\ (see, e.g., Ref. [H, p. 
751). The arbitrary constant in Eq. Q1Q) has been chosen 
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to satisfy the required asymptotic behavior V(z) — » (z + 
const) 2 /2 at z — > — oo. Notice that, as V > and T4 < 0, 
we should work with the negative branch of the product 
log function: rj < and W(rj) < 0. 

Integrating Eq. (10), we obtain the traveling wave so- 



lution in an implicit form: 

J(V) + z + C = 0, 

where 

U(V) 



J(V) 



i-C 



-c ' 



(ii) 



(12) 



U(V) = - In [— W (- exp(-l - V" 1 ))] and C is an ar- 
bitrary constant. 

To understand the asymptotic behavior of this solu- 
tion at z — > — oo and z — » +oo, we need to know the 
asymptotics of J(V). After some algebra we obtain 



j(Y) 



lnV+ Ax +0(V) , 
l 

3 

at V — > +0 and F — > +oo, respectively. Here 



(21/) 1 /2 + l\nV + A 2 +(D(V- 1 / 2 ) , 



(13) 



Ai = 



(1 - e-C) dC 
C(l-C-e-C) 



= -1.46074400... 



A 2 = -2 - (l/3)ln(2e) - A 3 , and 
l 

A 3 = 



1 



i _^_ e _ f + iL + A] ./, = -u.or,H,ns ! ,2... 



Using Eqs. (O) and (h_3!) we obtain: 



-C-Ax 



+ O (e _2z ) at z — > 4-oo (14) 



and 



U(z) = (l/2)(z + C + A 2 ) 2 

+ (2/3)(z + C + A 2 ) In (jz + C + A 2 \)/V~2) 



■0(\n 2 \z + C + A 2 |) at 



(15) 



The required asymptotic behavior V — > c exp(— z) at 
z — > -(-oo selects C = — Ai — lnc, so the traveling wave 
solution ( |lT| ) is now fully determined. After some rear- 
rangement, we rewrite the asymptotics dlj) and (|l5|) as 



V(«) = ce- z + 0{e~ 2z ) 



at 



z — > +oo . 



(16) 



and 



t/(z) = (1/2) (z - A) 2 + (2/3) (z - A) In \z - A| 

+ C(ln 2 |z|) at z^-oo, (17) 



where A = Ai - A 2 + (1/3) In 2 + In c. 

The leading term in Eq. (|l^) corresponds to a steady- 
state solution in the physical variables, while the sub- 
leading term is exponentially small with respect to the 
leading one. This essentially static ("glassy") behavior 
of the solution at large distances reflects effective diffu- 
sion choking at small temperatures. 

Now let us return to the bubble core description, Eq. 
(Q). Our basic assumption here (supported by the re- 
sults) is that, in the asymptotic stage r» 1, the terms 
Mo, mi, . . . depend on time only through the time depen- 
dences of a, of a and d, of a,d, respectively. The 
small parameter of this expansion is a/ a. In the zeroth 
approximation of this perturbation scheme, uq obeys Eq. 
(H) without the time derivative term. In the first approx- 
imation we obtain the following linear equation: 



L m (£) = ad cos 2 £ (1 + £ tan £) 



(18) 



where we have again used £ = x/a. 

The zero modes of the operator L are T(£) = cos 2 £ + 
£ cos£ sin£ and <&(£) = sin£ cos£. Looking for the 
general solution of Eq. (|l8| ) in the form of ui = 
Ci(0 T (0 + C 2(£)$(£) and defining a(r) by the con- 
dition u(0,t) — a 2 (r)/2, we arrive at 

= -(2/3) ad cos 2 £[(£tan£- 1) lncos£ 

+ 21n(2/e)£tan£ + tan£Im [Li 2 (-e 2i? )]] , (19) 

where Li 2 (x) = 53fcLi k~ 2 x k is the dilogarithm (see, e.g., 
Ref. p|, p. 743). In the vicinity of £ = tt/2 



ui = J ad | In ^ + ad O f |£| 3 In ||| 
3 V 



(20) 



where £ = £ — tt/2. Too close to £ = tt/2 the correction 
Mi and its derivatives become larger than the zero-order 
solution Mo and its corresponding derivatives, so the per- 
turbation procedure breaks down. Therefore, the bubble 
core solution [Eqs. @, (0) and (|n|)] should be matched 
with the traveling wave solution [Eq. (pr|)] in the region 
where \x — ira/2\ is small enough (so that the leading 
term of the asymptotics of Mo is much larger than the 
subleading terms) but, on the other hand, large enough 
(so that Mi is small compared to Mo). Working in this 
region and collecting the leading contributions from mo 
and Mi, we obtain after some rearrangement: 



= I & + ^d x In |£| + \ O (£ 4 ) + ^ O (\x\ 3 In \x\) 



C(d 2 ln 2 a) 0(\n\x\) 



(21) 



where x = x — (tt/2) a — (tt/3) a ln(e 2 a/4). 

Now we can perform the matching procedure. We re- 
quire that the z —* — oo asymptotics of the traveling wave 
solution, Eq. ([l7]), coincide with the asymptotics ( pl| ) of 
the bubble core solution. This yields 



3 



a(r) 



2 



2 1 T 
t In — 

3 4tt 



S + lncfc 2 + o(l) 



(22) 



where B = 1 + A x + A 3 = -0.514362926 . . ., o(l) denotes 
terms that vanish as r — ► oo and we have restored the 
fc-dependence jll|. We see that the leading term in a(r) 
is logarithmic in physical time t, and it coincides with 
Eq. (||). The subleading term behaves like lnr ~ lnlnt, 
while the sub-subleading term is constant. 

The matching region is determined by the require- 
ments that the subleading term in Eq. ( pi] ) is much 
less than the leading term, but much greater than the 
rest of terms. These yield an approximate condition 
In 2 a -C —x + | a <C a 2 / 3 , so that the matching region 
expands as t — > oo. 



for moderate Mach numbers, with those obtained with 
Eq. (P . This shows robustness of the reduced equation 
in the description of CCFs. 

In conclusion, we have shown that hot bubbles in gases 
expand logarithmically slowly in the process of cooling. 
By constructing an asymptotic solution, that matches 
a "quasi-similarity" inner solution and a "glassy" outer 
solution, we have been able to see how logarithmic cor- 
rections enter dynamic scaling. 

We are grateful to Y. Kurzweil for help with Fig. 1. 
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FIG. 1. The amplitude a versus the new time r = In £ found 
numerically (dashed line) and analytically (solid line). The 
parameters are described in the text. 

We compared the asymptotic solution with numerical 
simulations. Finding logarithmic corrections numerically 
usually requires going to very long times. Instead, we 
directly solved Eq. (^) in the new variables, which en- 
abled us to reach r ~ 20, that is t ~ 5 • 10 s . Eq. (g) 
was solved on the interval x £ (0, 9) subject to the Neu- 
mann boundary conditions. The initial condition was 
u(x, t = 0) = 2 cxp [-3 {x 2 + 0.1) 1 / 2 ] , so that c = 2 and 
k = 3. The system length was large enough for the solu- 
tion to enter the asymptotic regime before the expanding 
"core" reaches the boundary x = 9. Fig. 1 compares 

1/2 

a = (2u(0,t)) found numerically with the prediction 
ofEq. (p2|). At long times the agreement is excellent. We 
also verified other attributes of the asymptotic solution. 

In separate simulations that will be presented else- 
where Jl3| , evolution of the same initial condition was in- 
vestigated in the framework of the full gasdynamic equa- 
tions. Remarkably, the results essentially coincide, even 
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